Hydrodynamic relation in 2D Heisenberg antiferromagnet in a field 
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The spin-stiffness p 3 of a 2D Heisenberg antiferromagnet depends non-analytically on external 
magnetic field. We demonstrate that the hydrodynamic relation between p s , the uniform suscepti- 
bility x, an d the spin-wave velocity c is not violated by such a behavior because similar non-analytic 
terms from all three quantities mutually cancel out. In this work, explicit expressions for the field- 
dependent spin stiffness and for the magnon velocity of the 2D square lattice antiferromagnet are 
obtained by direct calculation to order 1/S and in the whole range of magnetic fields. 
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The effective description of spin- waves in the Heisen- 
berg and easy-plane antiferromagnets by a hydrody- 
namic theory goes back to the work by Halperin and 
Hohenbergii Such a description implies the following hy- 
drodynamic relation 



Ps 



= 1, 



(1) 



between the susceptibility x, the spin-wave velocity c, 
and the spin stiffness p s . The importance of an indepen- 
dent verification of such a relation using direct micro- 
scopic calculations has been recognized and received a 
significant attention in the past. 2 Corresponding calcula- 
tions confirming the validity of such a relation for the 2D 
square-lattice Heisenberg antiferromagnet (HAF) have 
been carried out in the early 1990s using the spin-wave 
theory^ to orders 1/S 2 and 1/S 3 . Numerical studies 2 ^ 
of the S — 1/2 case of the same model have also given 
a strong support of the relation (JTJ) . While initial inter- 
est in this problem was motivated by the large- J high- 
T c materials, more recently, synthesis of small- J quan- 
tum antiferromagnets^ has generated significant interest 
in the effects of external magnetic field in the properties 
of the HAFs, the regime that was previously unreachable. 

Uniform magnetic field lowers the full rotational sym- 
metry of the Heisenberg model to 0(2), making it equiv- 
alent to that of the the easy-plane antiferromagnets with 
the easy-plane of spin rotations perpendicular to the di- 
rection of the field. Note that the hydrodynamic consid- 
eration of Ref. U is also valid for the easy-plane antifer- 
romagnets. Thus, at the first glance, it seems natural to 
assume that the two hydrodynamic descriptions should 
connect continuously. However, the situation is far less 
trivial as several quantities were shown to exhibit a non- 
analytic behavior in small fields. In the earlier work, 
Ref. [A, field dependence of the ground-state energy and 
susceptibility was discussed for the non-linear er-modcl. 
The non-analytic field-dependent corrections have been 
found in the dimensions D < 3. A subsequent inde- 
pendent study, Ref. HI, obtained the same non-analytic 



behavior in small fields in the framework of the spin- 
wave theory. The recent work, Ref. 0, used a hybrid 
1/S-expansion— <x- model approach to demonstrate that 
the spin-wave velocity in a 2D antiferromagnet also has 
a non-analytic dependence on the field, c(H) — c(0) oc \H\ 
in the first 1/5 order. Recent studies of the combined ef- 
fects of the Dzyaloshinskii-Moriya and uniform magnetic 
field in the spectrum and the ground-state propertie d 10 ! 11 
of the 2D HAFs have also found non-analytic dependen- 
cies that are related to the ones discussed here. 

The origin of the non-analytic behavior can be traced 
to the field-induced gap in one of the Goldstone modes. 
External field creates the so-called uniform-precession 
mode, which corresponds to the precession of the field- 
induced magnetization around the field direction with the 
energy equal to H. When the field is small, the mode 
is almost gapless and contributes to the fluctuation cor- 
rections to various quantities. These fluctuations may, 
potentially, induce non-hydrodynamic corrections in the 
corresponding 1/S order of the theory. Thus, the validity 
of the relation ([TJ) in a field has to be verified. 

In the case of the square-lattice HAF, out of three con- 
stants needed for the hydrodynamic relation it is only 
the spin stiffness for which the presence or absence of the 
non-analytic terms in the field-dependence remains un- 
known. In this work we carry out direct analytical calcu- 
lations of p s to the necessary order in both 1/S and H to: 
(i) identify such non-analytic terms, and (ii) verify that 
the non-analytic behavior of all three quantities does not 
lead to the violation of the hydrodynamic relation (Q]) . In 
the course of such derivation, we also obtain an analytic 
expression for the spin stiffness to order 1/S and for all 
ranges of the field. In addition, the non-analytic behavior 
of the spin-wave velocity, previously obtained in Ref. 9 
by a hybrid 1/S — er-model approach, is confirmed within 
the framework of the spin-wave theory and the compact 
analytic expression for the velocity renormalization is ob- 
tained for an arbitrary value of the field. 

We would like to make a separate note on the recent 
work, Ref. that combines a thorough numerical in- 
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FIG. 1: (Color online) Left: angle of twist ip in the xo — yo 
plane of spins in one sublattice with respect to the other. 
Right: field-induced canting by the angle 9. 



vestigation of the static and dynamic properties of the 
5=1/2 square-lattice HAF in a field with the spin- wave 
analysis of the problem and provides a comprehensive 
comparison of the results. While in this work the spin 
stiffness is evaluated within the 1/5* spin- wave approxi- 
mation, it is done by numerical differentiation of the en- 
ergy with respect to the twist angle, and no non-analytic 
behavior of p s vs. H is discussed. Also, the hydrody- 
namic relation is used to provide a better estimate of 
the spin- wave velocity within the 1/ S approach, but the 
validity of it is not verified^ 

We consider the spin-S* HAF on the square lattice in an 
external field along the zq axis of the laboratory reference 
frame with the Hamiltonian given by 



(ij) « 



(2) 



where (ij) refer to the nearest-neighbor bonds. To study 
the spin stiffness, the Hamiltonian should be modified 
to introduce a twist angle between spins, rigidity to 
which should yields the stiffness directly. One of the 
prescriptions* is to twist spins in every second row by 
the fixed angle p. Another, intuitively more symmetric 
approach is to twist all the spins in one sublattice rela- 
tive to the other ^ In the latter method the twist energy 
is two times larger than in the former case because ev- 
ery spin has twice as many nearest neighbors that are 
twisted. For the Heisenberg model on a bi-partite lat- 
tice in zero field the direction of such a uniform twist is 
arbitrary. In the case of a non-zero external field, such 
a twist should be made in the plane perpendicular to 
the direction of the field, that is, in the xq — yo plane, 
see Fig. [1] Thus, using the sublattice twist with a small 
angle p <C 1, the modified Hamiltonian reads: 



H 



(ij) 



(3) 



where S^" = (Sf° , Sf°, 0) and we have omitted the terms 
that are linear in ip, as they either vanish or contribute to 
the p s -term only in the higher (1/S 2 ) order As such, 
the Hamiltonian © contains all the necessary terms to 
study both the classical limit of the model and the 1/S 
fluctuation corrections to it. 

To study quantum fluctuations around the classical 
spin configuration it is convenient to transform spins to 



"rotating" local reference frames in which the quantiza- 
tion axis z is along the classical spin direction4 ' 13 i 14 Mag- 
netic field cants spins toward its direction as is shown in 
Fig. [TJ Assuming that the spins lie in the x— z plane 
we perform transformation from the laboratory frame 
(xq, zq) into the rotating frame (x, z)&l& 

e iQri Sfcos9 , 



Sf 



= Sf sin ( 
= e lQr * Sf cos 9 + Sf sin 9 , Sf° = Sf 



(4) 



where Q = (tt,tt) is the ordering wave-vector, and cant- 
ing angle is as shown in Fig. [TJ The spin Hamiltonian ([3j) 
in the local coordinate system (fj| takes the form: 

H « J^2[SiSj -cos20 (S*S* + S*S*) 



(ij) 



(SfS] + S*S* sin 2 6 - S*S* cos 2 9) , (5) 



-H sinks'* 



where, again, the terms that are not contributing to the 
harmonic approximation are omitted. 

At the first glance, the spin stiffness can be defined 
from the averaging of the second line in Eq. ([5]) over 
the ground state. However, the situation is slightly 
more complex as the field-induced canting angle 9, which 
should be found from the minimization of the classical 
energy in ([5]), 8 also depends on the twist angled Per- 
forming such a minimization for {HJr* one obtains: 



sin ( 



(0) 



where the terms of higher order in the twist angle are 
truncated and the dimensionless variable h = H/(8JS), 
the field normalized to the saturation field H s = 8JS at 
which spins become fully aligned, is introduced. With 
the help of ([6]) one can eliminate 9 in ([5]) to obtain 



H 



Ti.ip=o + 7~Lp 



(7) 



where 7^=0 contains no twist angle and 7i Ps is given by: 

: J ^ \ V" 111 4- ^ '■ 2l 



7~i,, 



^2[(l + h 2 )S*S? + h 2 S*S 

m 



X QX 
j 



±sh 2 j2 s i 



(8) 



The subsequent treatment of the Hamiltonian 7i v= o 
involves standard bosonization of spin operators via the 
Holstein-Primakoff transformation to the first 1/5* order: 

Sf = 3- a\a t , » a^V2S , S+ = (S*") 1 , (9) 

which is followed by the Bogolyubov transformation. 8 
This yields the linear spin-wave theory Hamiltonian: 



n 



ip=0 



Egs 



4JS y^kQ^ak 



(10) 
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FIG. 2: (Color online) 1/S quantum correction to the spin 
stiffness ~Sp s / JS as a function of the field h. Inset: p a /JS 2 
vs h for S — 1/2, dashed line is pf/JS 2 vs h. 



inset presents p s for the spin-1/2 case. One of the inter- 
esting observations is that Sp s changes sign as a function 
of the field. It also exhibits a singular behavior in the 
derivative as h — *■ 1, similar to the one discussed before 
for the magnetization^ and is related to the logarithmi- 
cally vanishing scattering amplitude in the dilute 2D gas 
of bosons. The linear (non-analytic) field-dependence at 
small field is also clear from Fig. 

With the expressions (fT3"|) and (TTJJ at hand, one can 
now study the field dependence of p s at h — > 0. After 
some algebra one finds: 



jsj: 



i 



i 



,H=0 



+ 0(h 2 ). (15) 



It is easy to see that due to the field-induced gap cx H 
in the magnon spectrum the fluctuation terms like the 
one in (Tl5|) are yielding corrections oc \H\ in 2D. Some 
further algebra gives: 



see Ref. |8| for details. The dimensionless frequency is: 



v / (l + 7k)(l-(l-2/ l 2 ) 7k ) , 



(11) 



and 7k = (cosfc^ +cosfc y )/2. 

After diagonalization of Ti^o, the spin stiffness can be 
found as a coefficient in front of (p 2 in the twist part 7i p ^ 
of the Hamiltonian |8]) by averaging the spin operators 
over the spin-wave ground state and keeping terms to 
1/S order. Having in mind the extra factor of 2 in the 
sublattice twist approach, this finally yields: 



6E a 



(12) 



where p s = pf + Sp s with the classical and quantum 
contributions given by 

pf = JS 2 (l-h 2 ) , (13) 
d Ps = -JS[2n- A(l + h 2 ) + m(l - h 2 )] , 

where we use the following Hartree-Fock averages of the 
two-boson operator combinations: 



a\ai) 



{aiCij 



-E 

2 ^ 
q 



;(1 



h 2 ^ 



1 



to = (ajaj) 



It. 



7q 



(14) 



The above result (fl~3"|) gives p s to the order 1/S and 
for the fields anywhere between zero and the saturation 
value. At H = the expression for p s in (fl~3"|) and (fT4"]) co- 
incides with the known zero-field formula for the Heisen- 
berg modelfS^ The field-dependence of the quantum cor- 
rection to the spin stiffness Sp s is shown in Fig. O The 



Ps = JS 2 ( Z Ps + — \h\ )+0(h z ,l/S 2 ) 



„-H=0 



l + -Jy|/» 



(16) 



where p^ =0 contains zero-field, 1/S renormalization fac- 
tor Z Ps , Ref. [2|, and the last expression is obtained within 
the same 1/S* accuracy. 

For completeness, we also list the corresponding 1/S* 
expressions for the susceptibility. Magnetization of the 
square-lattice HAF is M = M c \ + 5M with the classical 
part and the quantum correction given by: 8 



M c i = Sh, SM 
Using x — dM/dH yields: 



-h(m 



1 



8y = — 
A 8J 



1 



i-^K 



X = Xd 

where I\ stands for the integral 



A - h 2 h) 



h =1^ 5 



(17) 



(18) 



(19) 



q "q 

At small fields, the same algebra as above gives: 

4 



X 



-(z x 



+ — \h\)+0{h 2 ,l/S 2 ) 



X 



H =°(l 



irS 



(20) 



Figure [3] shows the field-dependence of the quantum cor- 
rection to the susceptibility 5\- The inset presents \ f° r 
the case of S = 1/2. Similarly to Sp s , 8\ changes sign as 
a function of the field. As is discussed in Ref. d, \ has 
a singular logarithmic behavior at h — > 1 and the linear, 
non-analytic field-dependence at h — > is also clearly 
seen in Fig. [3J 
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FIG. 3: (Color online) 1/5 quantum correction to the suscep- 
tibility JS5x as a function of the field h. Inset: J\ vs h for 
5 = 1/2, dashed line is J\ci- 
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FIG. 4: (Color online) 1/5 quantum correction to the spin- 
wave velocity Ac/ J as a function of the field h. Inset: c/2JS 
vs h for 5 = 1/2, dashed line is cq/2JS vs ft. 



For the square lattice HAF in external magnetic field 
the energy of magnons to the first-order of the 1/5 ex- 
pansion is given byi£ 



e k = 4 J5w k + 6el 1] + 6e (2 



k ) 



(21) 



where Wk is defined in Eq. (fTTj) . correction 5e^ includes 
the Hartree-Fock and the canting angle renormalizations 

SejP = — \A-n + h 2 {A+m) [l - 2 7k (l-ft 2 )l (22) 
+ ~fl[n - A(l-ft 2 + 2ft 4 ) - rnft 2 (3-2ft 2 )] X 

( 2) 

and Se^ is the one-loop contributions from the three- 
magnon coupling: 



= -4Jft 2 (l-ft 2 ) 



(23) 



E 



Fi(k,q) 5 



F 2 (k,q) 2 



Wq+ Wk-q+Q— Wk Wk + Wq + Wk+q-Q 



Explicit expressions for ri(k, q) and r 2 (k, q) are given 
in Ref. [HI. After some algebra, the 1/5 correction to the 
spin-wave velocity can be written as: 



-co A(l-ft 2 +ft 4 )-n + mft 2 (2-ft 2 ) 
W ~ 5(1 - ft 2 ) 

h 2 2h 2 /I 2 (k) -4A 
~25 P 



(24) 



fc^O 

where the bare spin- wave velocity is Cq = 8J 2 5 2 (1 — h 2 ) 
and the k-dependent function in the last term is: 

t n \ on 7i«i/?2 + 72^1^2 , ocn 
J 2 (k) = 2(1 - ft )> 7i ? — , (25) 

' CJiW2(Wl + W2J 
q v 

with ai = (1 + 7i), A = [l - (1 - 2ft 2 )7 2 ] , and 1, 2 = 
q, q - k. 



The field-dependence of the quantum correction to the 
spin-wave velocity Ac, obtained by numerical evaluation 
of the integrals in Eq. ([24]) , is shown in Fig. [4] The inset 
presents the normalized magnon velocity for 5 = 1/2. 
The linear (non-analytic) field-dependence at small field 
is clearly visible in Fig. [4] Behavior of Ac at ft — > 1 is also 
singular, similarly to other quantities. It is interesting to 
note that the correction to the spin-wave velocity Ac is 
almost flat for 0.2 < ft < 0.9. 

It is easy to see that the first Hartree-Fock term in 
Eq. (|24p does not contribute to the anomalous non- 
analytic field dependence. After some more algebra, one 
can show that the same is true for the last term in (124p . 
The second term in Eq. (|24|) . on the other hand, yields 
at ft -> 0: 



Ac 

Co 



1 



1 



5 ^ VK) 3 «= ) 3 



JM 

irS 



This gives the same result as in Ref. S 



'■" ~ CH = ( 1 - \h\ 



(26) 



(27) 



Combining the expressions for the small-field expan- 
sion of all three quantities, p s , %, and c from Eqs. (| 16[) . 
(I2TJ)) . and [j2"7j). one can easily see that the hydrodynamic 
relation |T]) is obeyed as all the non-analytic terms explic- 
itly cancel each other in the leading order in ft. Moreover, 
such a verification of the relation {1} can be extended to 
an arbitrary field. Expanding the hydrodynamic relation 
(H|) to 1/5 order and observing that it is fulfilled at the 
classical level, XdCo/pf = 1, onc concludes that for the 
hydrodynamic relation to exist the quantum corrections 
from all three quantities must cancel each other at any 
field in each order of 1/5. For the 1/5 corrections this 
leads to 



Ps 



XclC 2 , 



^+2^-%)=0. ,,>S, 

Xcl c pf . 
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FIG. 5: (Color online) Cancellation of the 1/S quantum cor- 
rections in the hydrodynamic relation (|28p as a function of 
the field h (solid line). Dashed lines show contributions of 
individual terms in 



Numerical verification of the above relation is made using 
expressions (Hill) . ([T5]) . and (|2"4")1 and is presented by the 
solid line in Fig. [51 The dashed lines show contributions 
of individual terms in (f28|) . One can conclude that the 
cancellation takes place for all values of < h < 1. 

Having in mind the relation (|28[) . we can now obtain a 
much simpler expression for the spin-wave velocity rcnor- 



co 



Cf) 



-(- 



- n h 2 



(29) 



where 1\ is defined in Eq. (119)) . 

Altogether, we have confirmed the validity of the hy- 
drodynamic relation for the 2D Hciscnbcrg antiferromag- 
net in a uniform field. Despite the appearance of the non- 
analytic terms in the field-dependence of all key quan- 
tities due to quantum fluctuation involving small field- 
induced gap, they are not sufficient to violate such a 
relation. We have obtained expressions for the spin- 
stiffness p s and for the spin-wave velocity c for the square- 
lattice HAF, valid to the first-order in 1/S and for the 
whole range of magnetic fields. The non-analytic field- 
dependence of c, previously obtained by a hybrid 1/5- 
expansion— er-model approach, is verified using the more 
conventional spin- wave theory. 
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